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Abstract 

We use the method of the color effective Hamiltonian to study the structure 
of color singlet chain states in A'^c = 3 and in the large Nc limit. In order to 
obtain their total fraction when Nc is finite, we illustrate how to orthogonalize 
these non-orthogonal states. We give numerical results for the fraction of 
orthogonalized states in e"^e~ — > qqgg. With the help of a diagram technique, 
we derive their fraction up to 0{1/Nc) for the general multigluon process. For 
large Nc the singlet chain states correspond to well-defined color topologies. 
Therefore we may expect that the fraction of non-color-singlet-chain states is 
an estimate of the fraction of events where color reconnection is possible. In 
the case of soft gluon bremsstrahlung, we give an explicit form for the color 
effective Hamiltonian which leads to the dipole cascade formulation for parton 
showering in leading order in Nc. The next-to-leading order corrections are 
also given for e~^e~ qqgig2 and e+e~ qqgig2g3- 

(February 1, 2008) 

I. INTRODUCTION 

Hadronic processes in various high energy collisions are generally described in terms 
of two distinct phases: the perturbative phase and the non-perturbative hadronization 
one. The perturbative phase is well described by perturbative QCD (PQCD) while the 
hadronization phase cannot be described from first principle and can only be described by 
phenomenological models, e.g. the Lund string model or the cluster model 0. These 



*qwang@public. sdu.edu.cn 
^ gost a@ thep .lu.se 
■ffpcQsdu. edu.cn 
§xie@sdu.edu.cn 



1 



two phases are usually assumed to be well separated from each other. It is believed that the 
cross section for the hadronic process is fully determined by the perturbative phase, while in 
the hadronization phase a definite hadron state is chosen with total probability 1. In both 
phases, however, the large (the number of color) approximation is implied, which reduces 
the possible interference effects. A color charge of one parton is specifically connected to its 
anti-color in an accompanying parton, and with infinitely many colors the probability that 
two (or more) partons have the same color is zero So here enters the phenomenological 
color flow method (CFM) commonly used in the Lund model and the cluster model which 
implies assigning the color connection for the final parton system In these models, for 
every e^e^ —>■ qq + ng event, a neutral color flow is definitely determined and begins at the 
quark, connects each gluon one by one in a certain order, and ends at the antiquark. Each 
flow piece spanned between two partons is color-neutral and its hadronization is treated in 
a way similar to an independent qq singlet system. The present hadronization models work 
well, which shows that the CFM or the large Nc limit reflects some feature of the real world. 

Recently we proposed a strict PQCD method to study the color structure of a multiparton 
system. This method is called the method of color effective Hamiltonian which is constructed 
from the PQCD invariant amplitude. In this method a gluon is treated as an exact color 
octet, not a bi-color or nonet. The goal of this paper is to study the structure of the color 
singlet chain state (SCS or chain state) using this method for = 3 and in the large 
limit. 

The paper is organized as follows. In section |I|, we outline the effective Hamiltonian 
method. In section |T|, we define SCS, and introduce a diagram method which is very 
efficient in calculating the color inner product of any two color states. In section |V|, we 
calculate the fraction of SCS for Nc = 3. Since these states are not orthogonal to each 
other when A^^ is finite, we must find orthogonalized states to obtain a correct result. We 
introduce recipes for orthogonalization and use the orthogonalized states to calculate their 
fraction. As an example, we gives the numerical result for the fraction of orthogonalized 
SCS in the process of e~^e~ annihilation into quark, antiquark and two gluons. Section 
presents properties of SCS in the large limit. With the help of the diagram technique, we 
give their fraction up to 0{1/N^) for general multigluon process. In section [VT] , we give an 
explicit form for the momentum function D in He for the case of soft gluon bremsstrahlung. 
We show in the large Nc limit the consistency of He with the dipole cascade formulation for 
parton showering. We also discuss the next-to-leading order corrections for e^e~ qqgig2 
and e+e" qqgig2gz- 



II. COLOR EFFECTIVE HAMILTONIAN 

The QCD Lagrangian describes the SU{?>)c gauge interaction of gluon fields (a = 
1, ■ ■ ■ , 8) and quark fields q with three colors (i?, Y, B). By redefining Gell-Mann Matrices 
and gluon fields, we obtain eight new fields {Xf, Xf, Xf, Xf, Xf, Xf, X^, X^) which 
couple to color combinations {YR, RY, BY, YB, RB, BR, RR/YY, RR/YY/BB) re- 
spectively. Hence we write the QCD Lagrangian in a form with the quark-gluon interaction 
term showing clear color significance. This provides a strict formulation for calculating the 
fraction of color singlets for a multiparton system at the tree level from PQCD For the 
process e~^e~ — qq + ng, the essential part of the formulation is to exploit the color effective 
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Hamiltonian to compute the amplitude (/| He |0) for a certain color state |/). The color 
effective Hamiltonian He is found from the invariant amplitude M^^^"^ to be: 

where and ^^"l" are color operators for the quark and antiquark; (Q^)^ = is the nonet 

tensor operator; is the gluon's color octet operator; is a function of parton momenta 
and its dependence on the order of partons is marked by P which denotes a permutation 
of parton labels (1,2, The color effective Hamiltonian is another expression for 

the 5* matrix, and is therefore not Hermitian. One can verify its validity by using He to 
calculate the matrix element for the process e^e~ qq + ng, which returns the original 
|^aia2- an|2^ In Order to make He useful, we can define color states for the parton system 
qq + ng independently. A color state can be defined in the color space where each gluon 
subspace is either 8 or a little larger 3 (g) 3*. The effects of the unphysical " singlet-gluon" 
state, brought by enlarging the color space 8 to 3 (g) 3*(= 1 ® 8), can be eliminated by 
calculating the projection of the state on He- There are many ways of reducing the color 
space. Corresponding to each reduction recepe there is one set of orthogonal singlets spaces 
whose bases make up a complete and orthogonal set of color singlets. He has the property 
from unitarity, that the sum of the cross sections over all color singlets in a complete and 
orthogonal singlet set for the system qq + ng, is equal to the total cross section of e'^e" — > 
qq + ng at the tree level. 



III. SINGLET CHAIN STATES IN QQ + NG SYSTEM 



A chain state for a. qq + ng system is made up of a chain of (n + 1) pieces. Each piece 
is a singlet formed by the color charge of one parton and its anti-charge of the other one, 
and its color structure is the same as a qq system. A neutral color flow is also composed of 
pieces, each of which is formed by the color charge of one parton and its anti-charge of the 
other one, but it is a neutral color state not necessarily a color singlet. A singlet chain state 
can be decomposed into neutral color flow states. There are n\ chain states which connect q 
via n gluons in a specifled order to q, where the order of gluons is denoted by a permutation 
of (1, 2, • • • , n): P(l), -P(2), • • • , P{n). These states can be written as 

(I/,) , Z = 1, 2, • • • , n!) = {Ar-("+l)/2 |l,p(i)lp(i)p(2)lp(2)P(3) • • • lp(n)q)} (2) 

where A^~('*+^)/2 is a normahzation factor. Any two different chain states, \fi) and are 
not orthogonal to each other, i.e. (/j | /,) ^ 0. However, we can prove that any two different 
chain states are approximately orthogonal to each other to order As an example, we 

look at the inner product of two states as follows: 

l/l) = Ar-(-+l)/2 |l,ili2l23l3ll45 ' ' ' h,+l ' ' ' Inq) 
1/2) = iV-(-+l)/2 |l,ili3l32l24l45 " " " h,i+l ' ' ' Ing) 
C12 = (/l I /2) 

= A^^ '■"^^''(Iglll2l23l34l45 " " " " " " ^nq \ Iglll3l32l24l45 " " " " " " ^nq) 
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The gluon order in and I/2) is (12345 ■ ■ ■ ■ ■ - n) and (13245 • • • zz+1 • • ■ n) respectively. 
The only difference is that the position of (?2 and (73 is interchanged in two states. We write 
the explicit form of and I/2) as: 

/i) = A/'-("+^)/2^(ao, 6i)5(ai, 62)^(02, 63)^(03, 64)^(04, ^5) ■ ■ ■ 5(ai, 6^+1) ■ ■ ■ 5(a„, 60) 

V]/ Vl/^l\]> \]>^2x[/ vl/''3v]/ V|/64x]/ v[/fc5...vl/ Vl/''i+l . . . V|/ V]>''0 

flo f^-l C12 '2-3 '^4 f^-ri 

/2) = iV-("+i)/25(a[„ 6'3)5(a;„ 6'2)5(a'2, &'4)5«, ^'5) • • • • ■ ■ &o) 



where a,b,a',b' are the color and anticolor indices, and 0, 1, 2, 3, ■ ■ ■ , n denote q{q), gi, (72, 
93,- ■ -^Qn- Thus (/i I /2) is: 

(/i I /2) = A/'~^''+^H5(ao,6i)5(ai, 62)5(02, &3)'5(a3,&4)<5(a4,&5) ■ ■ -5(0^,6^+1) ■ ■ ■(5(a„,6o)] 
■[5(a'o, 6;)5(a;, 6'3)5(a'3, 6'2)5(a'2, &0<^(«4, ^s) " " " ^W, ^^+1) " " " ^^o)] (3) 

When n?i=o,i, --," ^i^ui ('''■u)^(pui is contracted with the content in the second square bracket, 
the above equation becomes 

(/i I /2) = A^^7(''+^)[5(ao,6i)5(ai, 62)^(^2, &3)'^(a3,&4)5(a4, 65) ■ ■ -5(0^,6^+1) ■ ■ ■(5(a„,6o)] 
■[5(ao, 61)5(01,63)5(03, 62)^(^2, &4)5(a4, 65) ■ ■ -5(0^,6^+1) ■ ■ ■5(a„,6o)] (4) 

We can use diagrams to visualize and simplify our calculation of inner products. Let us 
write the color and anticolor indices into two rows where the numbers in the first one are 
anticolor indices and those in the second are color ones. We draw a line between the number 
V in the first row and the number u in the second row if there is a 5( a^t, 6„) in (/, | fj). As a 
rule, we have (/j | fj) = Njr"'~^ where I is the number of closed paths. We can verify eq.(^ 
by drawing the corresponding diagram. As an example of (^), we show the I = 2 case in 
Fig.l. 

Generally, we can carry out the inner product of any two states by a diagram: 

{f^ I f)) 

= ^"+^^[5(ao,6p{i))5(ap(i),6p(2))5(ap(2),6p(3)) ■ ■ ■ 5(ap(i), 6p(i+i)) ■ ■■ 5(ap(„), 60)] 
•[5(ao, 6p'{i))5(ap/(i), 6p'(2))5(ap'(2), 6p/(3)) ■ ■ ■ 5(ap'(j), 6p/(j+i)) • • ■ 5(ap/(„), 69)] 

If there is a factor 5(ap(j), 6p(j+i))5(ap/(j), 6p/(j_|_i)) where P{i) = P'{j) and P{i + 1) = 
P'{j + 1), we get one factor Nc- n2 denotes the number of such identical 5-symbols in two 
square brackets of (|^). If we exclude the identical 5-symbols in the two brackets, the rest 
contains only different 5s. We denote the number of such 5s in one bracket as Ud- It can 
be demonstrated from the diagram that when rid is odd, these different 5s form one loop, 
and when rid is even they form two loops, rii denotes the number of loops formed by those 
different 5s and is 2/1 for even/odd n^. The maximum of ni + n2 — {n + 1) is —2, which 
occurs when ^2 reaches its maximum value n — 2. In this case the only difference between 
the two states is that only two gluons have their positions interchanged, just what we see in 
(§) and (^. Hence we see that any two different chain states are approximately orthogonal 
to each other to the order of 
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IV. ORTHOGONALIZATION AND FRACTION OF SINGLET CHAIN STATES 

AT FINITE Nc 



For finite Nc-, as in Nature, the chain states in (Q) are not orthogonal to each other. We 
cannot directly take the sum of each \{fi\H^\{))\' to give the total fraction for chain states, 
because it would make the contribution from the overlapped part of any two different states 
counted multiply. In this section we find a set of orthogonal states based on 

For the process e^e~ — >• qqgi ■■■(?«, there are n! chain states each connecting q 

with q through n gluons in a specific order. A set of orthogonal states {|//)} is related to 
{\fj)} by a linear transformation U}j, i.e. |//) = Ulj\fj), where Ulj is not a unitary matrix 
because of the non-orthogonality of the original states {\fj)}- U}j is not unique, so one 
can find many different ways of orthogonalization. Two different matrices and f/^ are 
associated with each other via a unitary matrix f/^^, which guarantees the conservation of 
probability. This means that for two sets of orthogonal states, \fl) and \ff)i which are both 
associated with the same set of non-orthogonal states the following identity holds: 

Due the rapidly growing complexity of orthogonalization for large numbers of emitted 
gluons, let us only discuss the two simplest cases e^e~ — > qqgig2 and e~^e~ qqgig2g?,- 
For e^e~ qqgig2, there are two chain states: 

(l/l) , 1/2)) = |l,lll2l2,) , ^ II.2I21I1,)) 

where and I/2) are not orthogonal to each other: (/i | = l/N^ = 1/9. A straightfor- 
ward way to construct orthogonal states is to linearly transform and I/2) into symmetric 
and anti-symmetric states respectively, i.e. 




, 1/2) 



3 3 



(6) 



where |/{) is the symmetric and I/2) the anti-symmetric state. Using {fi\Hc) = ^^-D — 

^^D^^ and (/2I i?c) = ~g^^^^ + ^73"^^^' obtain the sum of squared projections as 
follows: 

= + 3)(|Dif + iZ^^f ) + 2(Jft - 3)ReiD^' ■ D'^*) (7) 
^ 4.45(|Di2|' + |Z}2i|2) _ 3.i^epi2 . ^21*) 

The invariant amplitude for e+e" qq9i92 is: M^"' = {T''^T''^)abD^^ + {T''^T"^)abD^^ ■ 
The total cross section of e~^e~ qqgig2 at the tree level is then: 

atreeie'-e- qqg.g,) = J dQ {MT^Y .ox 

The functions D^^ and D^^ correspond to different kinematical distributions, which in general 
have a rather limited overlap. Therefore the kinematic interference term proportional to 
Re{D^^ ■ D"^^*) is suppressed. This kinematic interference term can be calculated in 2nd 
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order perturbation theory (we will come to this later). The result depends on the kinematical 
configuration, e. g. expressed by a y-cut for the definition of the 4-jet events. If we neglect 
kinematic interference terms in (0) and (^), we obtain: 

„/ + - - cr^cA cr(e+e- ^ qqgig2 ^ SCS) 

P{e^e qqgig2 SCS ) = — r ~ 83% (9) 

atree{e+e- qqgm) 

where we mention again that SCS is the abbreviation for singlet chain states. 

In order to evaluate the approximation of dropping kinematic interference terms, let 
us calculate the fraction P(e"*'e~ qqgig2 SCS) exactly. For a higher order process 
with more gluons produced, the calculation is too complicated. We know there are 8 lowest 
order Feynman diagrams in e~^e~ —>■ qqgig2, containing 2 diagrams with a tri-gluon vertex. 
In the calculation we choose the Feynman gauge, and we replace the polarization sum 
Sa=i,2 Ca(p)^a(p) with —g^^i where the sum is taken over two transverse polarizations A = 
1,2. However, —g^^ equals to the sum over all four polarizations, including two unphysical 
ones. To cancel unphysical polarization states and guarantee unitarity, we should introduce 
two ghost diagrams. Of course, we may work in the physical gauge, where there is no ghost, 
and directly use physical polarizations for the gluons, but the calculation in the Feynman 
gauge is much simpler; see Ref. for details. Thus, including two ghost diagrams, we 
have 10 diagrams altogether. The fact that ghost diagrams don't interfere with the 8 gluon 
diagrams make them easier to deal with. When we calculate the square of the amplitude 
M = Ml + M2 + ... + Mio, we know that a non-interference term |Mj|2 has a color factor 
16/3, while an interference term MjM* has a color factor —2/3. For a tri-gluon and ghost- 
ghost-gluon vertex, we make substitution T"T^ — T''T" = if°-^'^T'^ to get the functions D^^ 
and Z)2i. According to Eq.(^ and (|]), we have: 







1^21 


|')-3.1i?e(Di2.D2i*)] 











We use a Monte Carlo method to do the phase space integration and evaluate this ratio. 
The result is shown in Fig. 4 as a function of a cutoff i/cut = {Pi + Pj) ■ [Pi + Pj)/ s where Pi 
is the 4-momentum of parton i and ^/s is the center-of-mass energy of the e^e~ collision, 
which we set to 91GeV. In the figure, we see that the rate decreases slowly, from 0.72 to 
0.67, as Ucut varies from 10~^ to 10~2. These values are smaller than the rate 0.83 obtained 
by neglecting the kinematic interference contribution, but the difference is not large. This 
implies that kinematic interference terms are less important than the non-interference terms. 

For the qqgig2g3 system, there are 6 singlet chain states , i = 1, 2 ■ • • 6) which connect 
g to g via three ordered gluons: 123, 231, 312, 213, 132, 321, respectively. These states are 
not orthogonal to each other. The inner product (/j | fj) is 1 for i = j and 1/9 for 
i j. Our goal is to find 6 orthogonal states from them. The new states are denoted 
(l/i) 5 ^ = 1, 2 ■ ■ ■ 6). Making use of the symmetric and approximately orthogonal properties 
of {fi I fj), we find one orthogonal set as follows: 

I/O = (1 + a) I/,) + 6^:1/,) , /orz = l,---,6 (11) 

where a ~ 0.021, e ~ —0.037. According to arguments given for e^e~ qqgig2, we may 
find a different set of orthogonal states which is related to 1/^') via a unitary transformation, 
and we know that either set is equivalent. 
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Now we try to calculate the total fraction for chain states I//) by projecting them on 
\Hc). Since each state is orthogonal to any other, we can sum up all squared projections: 

ELi M I ^c)|' = 5.47{|D123|2 ^ |^231|2 ^ |^312|2 

+ |/}213|2^ |^132|2^ |^321|2| (12) 

+kinematic interference terms 

where kinematic interference terms refer to interference terms between two different D. Thus 
we obtain for chain states in e~^e^ — > qqgig2g3- 



a(e e 



scs)= ldnj2U\HM^ (13) 

i=i 



In order to estimate the fraction of chain states, we need to know the total cross section for 
e"'"e~ — >• qqgig2g3 at the tree level: 



atreeie^e- ^ qqgM) = I dn ■ M^"^"^ ■ (M^ 
= J dfl{j:p{T''PWT"p(^'>T''p(^))abD^){J2p,{T''p'wT''p'(^)T"p'(^))abD 



Here we resume the convention that a repetition of indices stands for summation. Expanding 
(0) gives: 



atree{e^e qqgig2g3) 

= Jdn- {i(|^[|Z)123|2 ^ |^231|2 ^ |^312|2 
+ |/}213|2^ |^132|2^ |^321|2] 

+kinematic interference terms} 



(15) 



where kinematic interference terms are suppressed by powers of l/Nc- Thus, to leading order 
in Nc, we can give an instant estimate for the fraction from (pTsD and (^) without carrying 
out phase space integrals: 

+ _ a{e-^e-^SCS) 5.47 
P(e+e ^ SCS) = — ^ = = ~ 77% 
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2 {2N^) 

In this section, we have discussed the orthogonalization for chain states and estimated 
their fraction for e"'"e~ — > qqgig2 and e"'"e~ — qqgig2g3 by neglecting kinematic interference 
and then keeping only the interference due to the finite number of colors, Nc- There are 
many ways of constructing orthogonal states from the original non-orthogonal ones. Different 
ways lead to different orthogonal states. They are, however, equivalent for calculations of 
probabilities. Normally we can make use of the fact that the original states are approximately 
orthogonal up to 0{1/N^). Thus one can find a set of orthogonal states, which are slightly 
different from the original ones, i.e. the transformation matrix is close to the unit matrix. 
The other straightforward orthogonalization recipe we give in this section, is to symmetrize 
and anti-symmetrize non-orthogonal states. We know that kinematic interference terms are 
all suppressed by 0{1/Nc) with respect to the non-interference terms. If we neglect all 
kinematic interference terms and then keep only the color interference brought by finite Nc, 
the total fraction is 83% for e~^e~ qqgig2 and 77% for e"'"e~ — > qqgig2g3- For the sake 
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of estimating the magnitude of the kinematic interference, we give the numerical result for 
the fraction P{e~^e~ (1(19192 SCS) with the kinematic interference taken into account. 
The result is shown in Fig. 4 as a function of t/cut- We see that the rate decreases slowly, from 
0.72 to 0.67, as ycut varies from 10~^ to 10~^. These values are smaller than the rate 0.83 
obtained by neglecting kinematic interference terms, but the difference is not large. This 
implies that kinematic interference terms are less important than non-interference ones, 
though not negligibly small. 



V. PROPERTIES OF SINGLET CHAIN STATES FOR QQ + NG IN THE LARGE 

Nc LIMIT 



In this section, we will study the properties of the chain states for qq + ng in the large 
Nc limit and obtain their fraction to 0{1/N^). 
The projection of a chain state |/) on \Hc) is: 

(/ I H,) = J2{l/V2rD^{f I Tr{QGp^,)Gp^,) ■ ■ •Gp(.))) 
p 

where is given by Without loss of generality, we choose 

|/)=iV-("+^)/2|l^lll2l23...1.,+l...l„^) 

The projection of any other chain state can be obtained by permuting gluon labels. For 
convenience, we will ignore the normalization factor N'^"^^^^^"^ and denote |/) as equiva- 
lent to |lgili2l23 • ■ • ■ ■ ■ Ing) in intermediate steps of the calculation. We will put the 
normalization factor back in the final results. 

The order of the n gluons in |/) is (1, 2, 3, ■ ■ ■ , n). In \Hc), there is also a term 
\Tr{QGiG2 ■ ■ -Gn)) with gluon labels in the same order. Let us first calculate 

(/ I Tr{QG,G2 ■ ■ ■ G,,)) = iV-("+^)/'(l,ili2l23 ■ ■ ■ • " • I TriQGiG2 ■ ■ ■ G„)) (16) 
We can expand \Tr{QGiG2 ■ ■ ■ Gn)) as follows: 



\Tr{QGiG2 ■ ■ ■ Gn)) = \Tr{QG[G'2 ■ ■ ■ G'J) 



(17) 



where (f i, f2, ■ ■ ■ , fn-fc) is the supplementary set to {ui,U2, ■ ■ ■ , Uk) in (1, 2, 3, ■ ■ ■ , n), and it 
satisfies Vi < V2 < ■ ■ ■ < Vn-k, i-e. the relative order of these n—k gluons in (G'^_^G'^^ ■ ■ ■ G'^^ J 
remains the same as in {QG1G2 ■ ■ ■ Gn)', J2{ul,■■■,u^:} sums over all decompositions of 
(12, ■ ■ ■ ,n) into (ui, M2, ■ ■ ■ , Uk) and {vi, V2, - ■ ■ , Vn-k)- Note that \Tr{QG'iG'2 ■ ■ ■ G'^)) is just 
|lgili2l23 ■ ■ ■ li,i+i ■ ■ ■ Ing)- Thus \Tr{QGiG2 ■ " 'Gn)) cau also be written as: 



\Tr{QGiG2 ■ ■ ■ Gn)) — |lglll2l23 ■ ■ ■ ' ' ' ^nq) 



l(Jl)llt)lt)2lt)2f3 ■ ■ ■ lfifi+1 ■ ■ ■ lfn-fc<? 



We find immediately 

(I,ili2l23 ■ ■ ■ l.,+i ■ ■ ■ 1„5 I Tr{QG\G'^ ■ ■ ■ G'J) = iV, 



(18) 



n+l 

c 
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Let us calculate 



(Iqlll2l23 ■ ■ ■ ■ ■ ■ Ing | lMlglll2l23 ' ' ' 1u-2,«--i1u--1,m+1 " ' ' ^nq) , 

This is one of the terms in (|1^) with k = 1. Therefore we denote it (/ | /c = 1 term) in 
short, We draw a diagram as shown in Fig. 2 where corresponding to each 1st, there is a 
line starting from position s in the lower row to position t in the upper row. Counting the 
number of closed paths, we obtain 



(Iglll2l23 ■ ■ ■ ■ ■ ■ Ing | Iulglll2l23 " " " ^u-2,u-l^u-l,u+l " " " ^nq) 

It is a little more complicated to calculate 

{f \ k = 2 term) 

(Iglll2l23 ' ' ' ' ' ' ^nq \ ^ui^U2^^qvi^viV2^V2V^ ' ' ' '^ViVi^i ' ' ' -^^n— fc?)) 



(19) 



(20) 



There are two cases: one where ui and U2 are adjacent to each other, the other where ui 
and U2 are not adjacent. To see more clearly, we look at the two cases separately, and we 
show the corresponding diagrams in Fig.3. For the case when ui and U2 are adjacent, we 
immediately have 

{f\k = 2 term) = A^^+i-s+i = n^-^ 
When Ui and U2 are not neighbors, we find 

{f\k = 2 term) = Ar;^+i-2x2+2 _ ^n-i 

Thus we obtain the same value in both cases. 
As a matter of fact, for a general expression 

(/ I k term) ^2^) 

= (191112123 ' ' ' lj,i+l ' ' ' Ing I (l'Uilu2 ' ' ' ^u-^.) (^^qvi^viV2^V2V3 ' ' ' "^ViVi^i ' ' ' I'Un-fc?)) 

where = 3, 4, 5, ■ ■ ■ , n, there is a unique value regardless of whether ui, U2,---,Uk or part of 
them are neighbors. Let us distinguish two cases with and without adjacent parton labels. 
If no label in -Ui, ■ ■ ■ ? Uk is adjacent to another, we have 

(/ I k term) = iv;+i-2fe+fc = iv;^-'=+i 

where n + 1 — 2A; is the number of double-line loops and k is the number of closed paths 
involving ui,U2, ■ ■ ■ ,Uk- If there are m labels, each of which is adjacent to at least one other, 
m can be grouped into / non-adjacent segments where labels belonging to the same segment 
are continuous, i.e. m = J2\^i m^ where m^ is the number of labels in the i-th. segment. 
Hence we get the result (/ | k term) = where e is 

e = [n + 1 - 2(A; - m) - (mi + 1) - (m2 + 1) (m; + 1)] . 

+ [k -m]+l = n~k + l ^ ^ 

Here the first term is the number of double-line loops, the second term is the contribution 
from k — m separated labels and the third term is from / continuous segments. Finally we 
have. 
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lnii\Tr{QGiG2---Gn)) 
= N^^' + ELi C^ii-lf^Nr'^' (23) 

where = ^^-^^i^^, denotes the number of ways to pick k out of n labels. 

Now we start calculating the general inner product (/ | Tr{QG p^i^G p{2) ■ ■ ■ Gp{n)))- First 
we focus on one of the simplest cases: 

{f\Tr{QG,G,G2G,G,---Gn)) (24) 

where two adjacent gluon labels 2 and 3 are interchanged compared to (|16|). Similar to (p!7D, 
we expand | Tr{QGiGzG2GiG^ ■ ■ ■ Gn)) as follows: 

|Tr(QGiG'3G'2G4G5---G„)) 

= \Tr{QG',G'^G'^G',G',---G'^)) (25) 



According to (^), the first term of (p^) is: 

(Iglll2l23 ■ ■ ■ ■ ■ ■ Ing | Igl I13I32 124I45 ' ' ' ' ' ' ^nq) 

Secondly we consider the following k = 1 terms: 

(-l/iVe)(lglll2l23 ■ ■ ■ ■■■lng\ I2(l,lll3l34 " " " h,i+l " " " ^ng)) (26) 

and 

( — 1/Ai"c)(lglll2l23 ■ ■ ■ h,i+l ■ ■ ■ Ing | I3(l(?lll2l24 " " " " " " ^nq)) (27) 

where gluon 2 and 3 are picked out as singlet I2 and I3 respectively. Note that the only 
difference between \Tr{QGiG3G2GiG5 ■ --Gn)) and \Tr{QGiG2G3G4G5 ■ --Gn)) is that the 



order of gluon 2 and 3 is interchanged. Thus we also encounter (|26|) and (|27|) in calculating 
(/ I Tr{QGiG2 ■ ■ ■ Gn)). From (0) we see that and (^ give the same value -N'^~\ 
For other terms we have 

(-l/iVj(/ I I„(l,ill3l32l24l45 " " " l.-2,.-ll«-l,«+l " " " In^)) = -N^'^ (28) 

where m 7^ 2,3. We see that ( ^8]) is suppressed by an additional factor compared to 

(p6|) and (^Tj) . It is easy to show that the contribution from terms 

(~^)'' |l«il«2 ■ ■ ■ l^fe) Tr{QG'^^G'^^ ■ ■ 



c 



with k > 1 are suppressed at least by l/N"^ as compared to (PS]). Thus, up to the highest 
order, we have: 

(/ I Tr{QG^G3G2G^G^ ■ ■ ■ G^)) = -N^' + O^N^-^) (29) 



10 



There are (n — 1) trace terms in |-ffc), which differ from |/) only in the order of two 
adjacent gluons. Their inner products with |/) are given by (p9|). There is a set of trace terms 
for which the labels of two non-adjacent gluons are interchanged compared to |/). We can 
in the same way prove that the projection of |/) on these trace terms gives N^^^ + 0{N^^^) , 



which is in the same magnitude as ( |29D but opposite in sign. These are next-to-leading order 
terms (~ compared to the leading term of the order N^^^ in (|23|) . Remaining terms 

have inner products which are suppressed by more powers of l/N^. 

Corresponding to a chain state |/), we can classify all trace terms in into three 
groups. One is the leading term where the gluon order is the same as that in |/). The 
second group are the next-to-leading terms, where the order of two gluons is interchanged 
relative to |/). This group can be further classified into two subgroups according to whether 
the order of two adjacent or non- adjacent gluons is interchanged respectively. The third 
group contains higher order terms, in which the order of the gluons differs even more from 
I/). The leading term is denoted as L{f), the next-to-leading terms as NL{f) with NLi{f) 
and NL2{f) for the adjacent and non-adjacent case respectively, and higher terms as H{f). 
Hence, for the state 

\f)=N-^^-^'^/^\lg,ll2h3---hi+l---lng) 

we find that, up to next-to-leading order, the projection on \Hc) is given by 
(/ I H,) = Epil/V2rDP{f I rr(gGp(i)Gp(2) ■ ■ •Gp(„))) 

^ (1/ V2)" ■ iV("+i)/2 . [(1 - ^)D^in - ^ Y.P^^,^u) + ^ T.P^NL,u) D"] ^^^^ 

We know from the previous section that there are n! singlet chain states which are denoted 
{\fi),i = 1,2, ■ ■ ■ ,n\}. They are not orthogonal to each other. The largest inner product 
of two states is l/N^- Suppose we find a set of orthogonal states {\f'i)} from {\fi)}. Up to 
next-to-leading order, we assume can be written in this form: 

+ E I/.) (31) 

where NL'{fi) refers to the set of chain states which contribute to 1//) in next-to-leading 
order; Ci and C2 are constants of order 1. One can verify that the set NL is included in 
NL'. The reason is that for a chain state the colors of the quark-antiquark pair play an 
equal role as those of gluons, while for a trace term in \Hc) it emerges as a color nonet 
which is different from the color octets of gluons. To see it more clearly, we take as an 
example two permutations (01234) and (02341). Obviously for chain states (01234|02341) = 
(l/iV2)(01234|01234), hence the orders (01234) and (02341) belong to the same NL' set. 
Here we use a simplified notation for the chain state, e.g. |01234) = |loili2l23l34l4o)- It is 
easy to understand this because we can write the second state in the form 1 10234), which is 
different from the state |01234) in two labels and 1. We can verify that (01234|Tr(02341)) 
is suppressed by 0{1/N^) relative to (01234|Tr(01234)), i.e. |Tr(02341)) does not belong 
to A^L(01234). Here we use a simplified notation for the chain state, e.g. Tr (01234) = 
Tr{QGiG2G^Gi). Another different feature of NL' compared to NL is that the inner 
product of any two different states, where only two gluon labels are interchanged, is always 
suppressed by 1/iV^ relative to the inner product of themselves, regardless of whether they 
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belong to the adjacent or non- adjacent case. Hence NV can be written in this form: NV = 
NLi + NL2 + WL. Up to 0(1 /N^) we have 

(// I H,) = (1 + I H,) - ^ T^jeNLifM I He) 

+ ^ EpeNL^in) D"-^ E.eNLif.) D^^f^n (32) 
= 2-"/2iV("+^)/2[(l + ^)D^(f^) - ^ EpeNL.if.) 



and the projection square has the following form: 

|(/; I = 2-iVrM (1 + - ^ Ep.^l,(/.) i?e(D^(/>) ■ D^*) 

where we only keep the next-to-leading order. The total sum is 



(33) 



Er=ii(//i^e)p 

= 2-N-^\ (1 + ^) Ep - ^ E{P,P'}eiVL. i?e(D^ ■ D^'*) (34) 

E{P,P'KiVL. i?e(D^ ■ D^'*) - ^ Ep,ivL(/o ^e(D^(/') ■ D^*) ] 

where {P, P'} G A^Li means that the orders P and P' are different in only two neighboring 
gluon labels, {P, P'} G A^L2 means that they are different in two non-adjacent gluon labels, 
and {P, P'} G iVL means {P, P'} G A^L' but not included in NL^ + A^Lg. 

Now we start calculating crtr.ee(e^e~ — > gg + ng) in the large A^c limit. Recall that the 
ordinary matrix element is given by 



The total cross section is then: 



Ea,M.,..,...,«„/rff^|M:r"-"l' ^ (35) 



= Ea,M„a,,.,a„ /rff^ ^p(T-^(i)T°^(2) . . .T"^(n)),,Z}^ 
In evaluating (p5|), we mainly encounter two types of traces of Gell-Mann matrices 

']^r(T°'^T"'^ ■ ■ ■ X'""T"" ■ ■ ■ ^"'"a^"'"!^ 

= N,-C^ = N,{^r (36) 
^ (1/2)"- Ar^("+i)(l- ^) 



and 

TriA.T-A^T'-As) = '^Tr{A,)Tr{A^As) - ^Tr{A,A,A, 

lTr{A,A,) - ^ 



Tr(Air")rr(A2T-) = \Tr{A^A2) - 2^rr(Ai)Tr(A2) ^""'^ 



where Ai,A2 and ^3 are chains of products of Gell-Mann matrices. The result in (BBI) is just 



the color factor associated with the terms 



D 



p 



2 



while ( P?!) gives the color factor for the 
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kinematic interference terms D^-D^' , where P and P' denote two different permutations. In 
the large Nc hmit the factor (|36|) is the leading one. The next-to-leading contribution comes 
from ■ terms where P and P' are different in only two gluon labels. Also the next- 
to- leading contributions can be classified into two different cases, one is that {P, P'} G NLi, 
i.e. P and P' are different in two neighboring gluon labels, the other is that {P, P'} G NL2, 
i.e. they are different in two non-adjacent gluon labels. For {P, P'} G NLi, the color factor 
is ■ ■ C^-^ ^ -j^N^~\ while for {P, P'} G NL2, the color factor is j^N^'^ in the 
large Nc limit. After keeping terms up to next-to- leading order we have 

atree{e^e- qq + ng) 

= (1/2)- ■ Ar("+i) . / dn[ (1 -^^) Ep ID^' - ^ E{P,P'}eNL, ReiD^ ■ D^'*) (38) 
+ ^ T,{p,p'}eNL2 Re{D^ ■ *) ] 

We write the fraction of singlet chain states as: 

P(e+e- ^qq + ng SCS) . . 

= (/ dn Eti I {fl I Hc)\')/atree{e+e- qq + ng) ^^^^ 



where YJi=i \ {fi I -f^c)P is given by P^). According to Eq.(P^, ^) and we reach our final 
result up to next-to-leading order 

P(e+e- ^ qq + ng SCS) 

= 1 + ^[(2Ci -n)- (4C2 + 2)Ti - (4C2 - 2)T2 - 4C2T3] ^^^> 
where Ti and T2 are defined by: 



^ _ / E{P,P^}g^L. Pe(D^-D^'-) 

^ /^^^EpI^^I^ 

^ ^ /rf^^E{P,p.KlvIi?e(D^■D^'-) 

' /c^^^EpI^^I' 

The result in Eq.(^D|) gives the fraction of events, for which the colors of the gluons cor- 
respond to a single chain from the quark to the antiquark. We expect that these states 
hadronize producing a corresponding chain of hadrons. The remaining events correspond 
to more complicated color structures, where one gluon can be connected to more than two 
other gluons, as indicated in Fig. 2 and 3. Some dynamical feature of the confining mech- 
anism may imply that also these states result in string-like hadronic final states, but it is 
also conceivable that these parton states can produce more complex hadron configurations. 



VI. He IN SOFT GLUON BREMSSTRAHLUNG AND ITS RELATION TO 

DIPOLE CASCADE MODEL 

We have not yet touched the momentum function D in He so far. The number and com- 
plexity of Feynman diagrams increases drastically with growing number of emitted gluons 
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in normal situations. However, in the case of soft gluon bremsstrahlung, the D function 
has a simple and recursive structure. The method we use to derive the D function or 
in this section is called soft gluon insertion technique 0. By recursively adding a softer 
(with lower energy) gluon in multigluon emissions, the distribution for each new emission is 
approximately determined by an eikonal current stemming from all the harder gluons. The 
result factorizes between the emissions and is equivalent to classical bremsstrahlung under 
certain angular ordering conditions. 

Assume that in e~^e~ annihilation, e^e~ —>■ qqgig2 ■ ■ ■ Qn ■, n gluons are all soft ones and 
their energies/momenta are strongly ordered: 

~ Ep. > Efc, > Efc, ■ ■ ■ > (42) 

where p, p', ki,k2, - ■ ■ , kn are 4-momenta of q, q, Qi, 92, ■ ■ ■ , Qn respectively. When the hardest 
gluon Qi is emitted, it has two legs to attach to, one is the quark's momentum and the other 
is the anti-quark's. They give rise to the amplitude: 

~ 9se,M^, = gsi^ - -^)T2e,, = g,J^^{k,;p,p')T:,^e,, (43) 

P ■ Ki P ■ rCl 

where e denotes the gluon's polarization 4- vector and Qs is the strong coupling constant. 

When the second hardest gluon g2 is emitted, it has three legs to attach to: q, q and gi. 
The amplitude is: 

ab ~ 5's^Mi^M2-'afeaia2 

= gs^fii£fj.2'J^^iki;p,p')[^(T°'^T°'^)ab - fr^(T°'^T°'^)ab + fc^^Aoi/j^ffe] ^4) 
= g^,e,,e^,[J^^-{ki;p,p')J>^^{k2;p,k,){T-^T'-%, 



+ J'^i {ki;p, p') J''^ {k2;ki, p') (T"! T 



where T^^^^ = ifaia2^ is the generator of the adjoint representation of SU{3), and the rule 
ifc.ia2f}T'^ = T^'^T''^ - T'^^T''^ has been used. 

In the case of n-gluon emission, we can prove in the same way that the amplitude and 
its corresponding \Hc) can be written in this form: 

~ g^e^.e^, ■ ■ ■ £^JEp(T"^'^)T"^<^' ■ ■ -T-^c^Oab 

•J'^i(fci;p,p')</^'(fc2; k2h, k2e) ■ ■ ■ J^-ikn] knh, ^ne)] (45) 

\Hc) ~ fe,£^,---£^„Ep(l/v^)lTr(gGp(i)Gp(2)---Gp(„))) 

■J^''{k^;p,p')J>'^{k2] k2h, k2e) ■ ■ ■ J^"ik^; knh, Ke) 

where subscripts ih and ie (where i = 2, ■ ■ ■ ,n) are determined by the following procedure: 
in the sequence (0, -P(l), -P(2), • ■ ■ , P{n), 0) (where we imply = g at the head and = q 
at the end), find the position of i, take away all greater numbers in the sequence, the left- 
nearest neighbor to i is ih and its right-nearest neighbor is ie. Having this complete form of 
\Hc), we can calculate the fraction of any color state by projecting the state onto it. Here 
we are only interested in what happens to SCS in the large limit. We consider a chain 
state l/p) which corresponds to a specific order of gluons: (-P(l), -P(2), ■ ■ ■ , P{n)) where P 
denotes the permutation of {1,2, ■ ■ ■ ,n) . According to the former section, in the leading 
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order, the inner product (/p | only picks up the term with the same order of gluons in 
\Hc): 

{fp I He) ^ gs^f^^e^^ ■ ■■€^,J^^{kuP,p')Jf'^{k2] k2h, k2e) ■ ■ ■ J^"{kn] knh, kne) ..^x 

•(l/x/2)-(/p I Tr(gGp(i)Gp(2) ■ • • Gp(„))) ^ 
According to (0), the projection square is: 

I {fp I H,) |2= N:^^gl^{p,p\,{k2h, k2e)k, ■ ■ ■ {k^h, fc„e)fc„ (47) 

where the antenna term is defined by: 

{k-pi){k-p2) 

Hence we see that each gluon gi is associated with two harder gluons Qih and Qie which are 
nearest to its position in the sequence. 
As an example, we look at the case: 

If) = iV-("+l)/2 |l^ill2l23 ■ ■ ■ ln+1 ■ ■ ■ Ing) (49) 

with gluons' order (1, 2, ■ ■ ■ , n). The projection square is: 

I (/ I He) \'= K+'g'^-{p,p%,{k„p\,{k2,p%, ■ ■ ■ (A;„_i,p')fc„ (50) 
The cross section of e~^e^ qq + ng ^ \f) is: 



drr(p+p- -->\f\\— Arn+l 2n d^p (fp' yin _d 



(2-k)'^2E (2tt)'-^2E' '■'■1=1 {2TT)'J'2Ei C5]_~\ 

{p,p')kAki.p')k,{k2,p')k, ■ ■ ■ (fc„-i,p')fe.(27r)M(^H^e+ + Pe--[p + p' + i: h]) 

In the approximation of soft gluon bremsstrahlung, the above cross section can be written 
as: 



da{eU^ - I/)) = da{e^e~ - qq) YIU f^s^^j^f 
= cia(e+e--,g)nr^if«.ff^f^ 



{ki-ki-i){ki-p') {52) 



where the strong coupling constant = and we use the notation ko = p for convenience. 
Define two new variables for the gluon gi which are called the generalized rapidity and the 
transverse momentum: 

yi9^-lq - 9^) = IH'^) , k?r{g,.,q ^ g,) = ^^'■^^g'''^') (53) 

We can rewrite Eq. ( [^ ) in terms of y and pt' 

n 

da{e^e' |/)) = d(r{e^e' qq) J| da{gi^iq gi) (54) 

i=l 

where we define the emission rate for gi'. 
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da[gi_iq gi) = —asdy[g,_iq gij-pr, = T (^5) 

27r kT[gi^iq gi) 

Eq. ( p^ is the result for one chain state with a specific order. If we take into account all n\ 
states, Eg . (|5^) becomes: 



da{e+e- E"=i \fj)) = da{e^e- qq)da{qq -> gi) 

■ [dcr{qgi -> g2)da{qg2giq ^ gs-'-gn) (56) 
+dcr{giq g2)da{qgig2q ^ 5-3 ■ ■ ■ 9n)] 

Eq.(p5|, |56| ) is just the dipole radiation formula in the Lund Dipole Cascade Model In 
a more general situation the natural ordering variable is transverse momentum, and it is 
also possible to go beyond the eikonal approximationso that the Altarelli-Parisi splitting 
functions are properly reproduced for coUinear emissions. In ref. it is also demonstrated 
that the dipole formulation reproduces the angular ordering constraint due to soft gluon 
interference [|10| . This restricts interference effects and implies that gluon emission is a local 



process, in which only a limited set of gluons contribute to the emisson of a softer gluon 
(where soft should mean in the rest frame of the parent set). 

We have derived the Lund dipole cascade model from the He approach in the leading 
order. Now we begin to discuss the next-to-leading order. The general case in this order is 
rather complicated, and we therefore first discuss the simple cases with two and three gluons. 
The local character of the dipole cascade emission implies, however, that these results may 
be relevant also for a more general situation. 

For n = 2, i.e. when there are two gluons in the final state, we have two singlet 
chain states |/i,2) = {|012), |021)} (the normalization factor is A^~^/^). We have NL' = 
{(012), (021)}, A^Li = NL = NL' and NL2 = NT={}. There are two types of inner prod- 
ucts between a chain state and a trace state in \Hc): (012|Tr(012)) = N^—2Nc+^ ~ and 
(012|Tr(021)) = ~Nc+^ ~ —N^. Note that we have used a simplified notation for the chain 
state and trace state, e.g. |01234) = II01I12I23I34I40) and Tr(01234) = Tr{QGiG2G3G4) 
etc.. If we use orthogonalized chain states as shown in Eq. (|3TD , we have: 



1(012)') = (l + f^)|012)-^|021) 

1(021)') = (l + g)|021)-f 1012) ^"'^ 

From Eq.(0) and taking into account the normalization factor for chain states, we obtain: 
|((012)'|i7,)P + |((021)'|iJ,)|2 

= lN^[ (1 + 2^)(|D12|2 + |^21|2) _ 4(^^gpl2 . ^21*) ] 



iV3^f{(l + 2^)(0,0)i[ (1,0)2 + (0,1)2] 

^(^(o,o)i[ (0,0)2 -(0,1)2 -(1,0)2]} 



(58) 



where we have used a simplified notation for the antenna term, for example, (0, 1)2 = {p, ki)k^ 
and (1, 0)2 = {ki,p')k2 etc.. We see that the magnitude of (1, 0)2 and (0, 1)2 in leading order 
obtains a correction term which is of order There is also a negative term (0, 0)2 from 

the interference which is absent in leading order. This term corresponds to the emission rate 
of gluon 2 from the dipole qq. 
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For n = 3, there are six chain states: 



{!/,), ^ = 1 ■ ■ ■ 6} = {(0123), (0132), (0213), (0312), (0231), (0321)} (59) 

Since {fi\fj) = ^/N^ for any i ^ j, all these states belong to the same NL'. NLi and NL2 
are given by: 

A^Li,2(0123) = {(0213), (0132)}, {(0321)} 

A^Li'2(0132) = {(0312), (0123)}, {(0231)} 

iVLi;2(0213) = {(0123), (0231)}, {(0312)} 

ArLi,2(0312) = {(0132), (0321)}, {(0213)} ^ ^ 

ArLi^2(0231) = {(0321), (0213)}, {(0132)} 

ArLi'2(0321) = {(0231), (0312)}, {(0123)} 



and NL is obtained by NL' — NLi — NL2. We can write the above equation in another 
form. According to Eq.(^), the following pairs of orders for chain states belong to NLi, 
NL2 and NL, respectively: 



NLi ■ {(0123), (0213)}, {(0123), (0132)}, {(0132), (0312)} 
{(0213), (0231)}, {(0312), (0321)}, {(0231), (0321)} 



(61) 

NL2: {(0123), (0321)}, {(0132), (0231)}, {(0213), (0312)} (62) 



NL : {(0123), (0312)}, {(0123), (0231)}, {(0312), (0231)} 
{(0321), (0132)}, {(0321), (0213)}, {(0132), (0213)} 

According Eg . (|6T| - |63D , the leading contribution involves 



(63) 



Y^aii p l^^'P = 23(?6(0, 0)1 { (1, 0)2 [ (0, 1)3 + (1, 2)3 + (2, 0)3 
+ (0,1)2 [(0,2)3 +(2, 1)3 + (1,0)3 ] } 



(64) 



The next-to-leading contribution consists of three parts which are from NLi, NL2 and NL 
respectively. The NLi contribution corresponds to the following sums: 

Tp^^L,Re{DPDP'*) = 2gl (O^O)i { -/i(0, 1, 0|2)/i(_l, 0, 2|3) 

-/i(0, 1, 0|2)/i(M, 2|3) + 2(1, 0)2 /i(l, 2, 0|3) + 2(1, 0)2 /i(0, 1, 2|3) (65) 

+2(0, 1)2 /i(2, 1, 0|3) + 2(0, 1)2 /i(0, 2, 1|3) } 

where we the interference pattern Ii is defined by: 



h{gi, 92, g2.\gi) = {gi, 93) - {91,92)94 - {92,93)94 

h{93, 92, 9i\9i) = h{9i, 92, 93\9a) 
where 9i is the label for gluon i. The A^L2 contribution corresponds to 

Y.P^^L,R_e{D''DP'*) =2gl (0,0)i /i(0,l,0|2) 
■[/i(0,2,0|3) + /i(0,l,0|3)-2(l,2)3] 



(66) 



(67) 
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The NL contribution corresponds to 



= 2 (O^O)i { 2 (1, 0)2 /2(0, 1, 2, 0|3) + 2 (0, 1)2 h{Q. 2, 1, 0|3) _ (68) 
-Ji(0, 1, 0|2)[ 2, 0|3) + /i(0, 1, 2|3) + /i(0, 2, 1|3) + /i(2, 1, 0|3) ] } 

where the interference pattern I2 is defined by: 

^2(5^1, 92, 93, 9^95) = (91, 9i)g^ + {92, 9z)g^ " {9l,93)g^, " {92, 9i)g^ (69) 

Substituting the leading and next-to-leading contribution given in Eg. ( |64| - |6^ ) into Eg. (^) , 
we can obtain the cross section for chain states up to next-to-leading order. 

In a strongly ordered gluon cascade the emission amplitude has an eikonal form. In 
this case, we have shown that in leading order in N^, the cross section of e^e~ — >■ qq + ng 
can be decomposed into n\ independent or incoherent parts and each part represents the 
contribution from the chain state with a specific order of gluons. The cross section of a 
singlet chain state exactly corresponds to the emission of gluons by a specific seguence of 
dipoles. The seguence is determined by the gluon order of the chain state. The softest gluon 
is emitted independently and in an egual emission probability by all possible dipoles each of 
which is stretched by two adjacent harder gluons, and each dipole corresponds to a specific 
chain state. Hence to sum over all possible chain states is the same thing as to sum over 
all possible dipole seguences. In next-to-leading order, the emission probability aguires a 
0{1/N^) correction and may not be necessarily egual for different neighbor-gluon dipoles. 
Furthermore there appear dipoles stretched by non-adjacent gluons which are separated by 
only one extra gluon. The next-to-leading corrections arise from interferences of two D- 
functions with their gluon orders in the NLi^2 oi NL set. For two D-functions belonging 
to A^Li 2, only one interference pattern Ii is relevant. The other interference pattern I2 is 
associated with D-functions belonging to NL. Ii is related to three consecutive gluons and 
I2 to four gluons which are not necessarily consecutive. In summary, each next-to-leading 
correction brought by the interference of a pair of D-functions, say, and , can be 
regarded as a small perturbation to the dipole seguence which corresponds to |-D^P or p 
by having dipoles formed by non-adjacent gluons while keeping the rest of the dipole seguence 
same as that corresponding to |-D^P or \D^ p. Compared to higher-order corrections, the 
next-to-leading correction causes the least perturbation to the dipole seguence of the leading 
order. 

VII. SUMMARY AND CONCLUSION 

The phenomenological color fiow picture, commonly used in the Lund model and the 
cluster model, is to assign the color connection of a final parton system. In these models, for 
an e+e~ ^ qq + ng event, the neutral color flow is definitely determined and begins at the 
guark, connects each gluon one by one in a certain order, and ends at the antiguark. Each 
fiow piece spanned between two partons is color-neutral and its hadronization is treated in 
a way similar to a qq singlet system. The present hadronization models work successfully, 
which shows that this picture is a good approximation to the real world. In this paper we 
use the method of color Hamiltonian, a strict formulation developed from PQCD, to study 
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the structure of chain states in e^e~ qq + ng for finite Nc = 3 and in the large Nc limit. 
For large these states correspond to well-defined color topologies. They just correspond 
to the phenomenological neutral color flow. Therefore we may expect that the fraction of 
the non-chain state is an estimate of the fraction of events, where color reconnection is 
possible. It is also conceivable that color structures, where one gluon is connected to more 
than two other gluons, may result in more complicated hadron configurations. (The problem 
of finding experimental observables for these final states is, however, beyond the scope of 
this paper.) 

There are n\ chain states, each of which corresponds to a specific order of n gluons. When 
Nc is 3, as in the real world, any two different chain states are not orthogonal to each other. 
To derive the total fraction for chain states we must orthogonalize them. We give two types 
of recipes of orthogonalization: one is to symmetrize the original chain states, the other is to 
find the transformation matrix which slightly differs from the unit matrix, by exploiting the 
fact that every two different chain states are approximately orthogonal up to order As 
an example, we give the numerical result for the rate of chain states in e'^e" — > qqgig2- The 
result is shown as a function of the cutoff ycut- The rate decreases slowly, from 0.72 to 0.67, 
as Vcut variates from 10~^ to 10~^. These values are smaller than the rate 0.83 obtained 
by neglecting the kinematic interference contribution. The difference is not large, which 
implies that kinematic interference terms are less important than non-interference ones, but 
not negligibly small. Therefore we may expect that the fraction of non-chain states is an 
estimate of the fraction of events where color reconnection is possible. 

Similar to n\ singlet chain states, there are also n\ terms in the color Hamiltonian, where 
each term corresponds to an order between n gluons. Up to 0{-^) we prove, with the help 
of a diagram technique in section 111, that there exists a one-to-one correspondence between 
a chain state and the term in with the same order of gluon labels. This means that 
when computing the fraction of a chain state |/) with a specific order of gluon connections, 
up to 0{1/Nc) we only need to consider the contribution from the term in He where the 
gluon labels are in the same order, and those terms for which the order of gluons is most 
close to it. 

Finally we give the explicit form for the D function and He in a special case, the case 
of soft gluon bremsstrahlung. In soft gluon bremsstrahlung with gluons strongly ordered 
in energy or transverse momentum, the emission amplitude has an eikonal form. In this 
case, we have shown that in leading order in N^^ the cross section of e^e~ — > qq + ng 
can be decomposed into n\ independent or incoherent parts and each part represents the 
contribution from the chain state with a specific order of gluons. We also give the next-to- 
leading order corrections for e'^e" — > qqgig2 and e'^e" — > qqgig2g3- The corrections arise 
from interferences of two D-functions with their gluon orders in the NLi^2 or NL set. Each 
next-to-leading correction brought by the interference of a pair of D-functions, say, and 

, can be regarded as a small perturbation to the dipole sequence which corresponds to 
|Z}^P or p by having dipoles formed by non-adjacent gluons, while keeping the rest 
of the dipole sequence same as that corresponding to or \D^'\'^. The next-to-leading 

correction causes the least perturbation to the dipole sequence of the leading order compared 
to higher-order corrections. 
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FIGURES 




FIG. 2. The diagram for (I51I12I23 • • • Iv+i • • • Ing | l«lgil23 • • • lw-2,«-il«-i,w+i • • • Ing)- 
There are (n + 1 — 2) double- hne loops plus one loop connecting color u — 1, u and anticolor 
u, u + 1. 




FIG. 3. The diagram for (lgili2 • • • • • • Ing I l«ilw2(l«wl 

and U2 are are neighbors, (b) ui and U2 are not neighbors. 
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FIG. 4. The numerical result for the fraction of chain states in e+e annihilation into two 
gluons. 
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